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ABSTRACT 

Observations of low-mass stars reveal a variety of magnetic field topologies ranging 
from large-scale, axial dipoles to more complex magnetic fields. At the same time, 
three-dimensional spherical simulations of convectively driven dynamos reproduce a 
similar diversity, which is commonly obtained either with Boussinesq models or with 
more realistic models based on the anelastic approximation, which take into account 
the variation of the density with depth throughout the convection zone. Neverthe¬ 
less, a conclusion from different anelastic studies is that dipolar solutions seem more 
difficult to obtain as soon as substantial stratifications are considered. In this paper, 
we aim at clarifying this point by investigating in more detail the influence of the 
density stratification on dipolar dynamos. To that end, we rely on a systematic pa¬ 
rameter study that allows us to clearly follow the evolution of the stability domain of 
the dipolar branch as the density stratification is increased. The impact of the den¬ 
sity stratification both on the dynamo onset and the dipole collapse is discussed and 
compared to previous Boussinesq results. Furthermore, our study indicates that the 
loss of the dipolar branch does not ensue from a specific modification of the dynamo 
mechanisms related to the background stratification, but could instead result from 
a bias as our observations naturally favour a certain domain in the parameter space 
characterized by moderate values of the Ekman number, owing to current computa¬ 
tional limitations. Moreover, we also show that the critical magnetic Reynolds number 
of the dipolar branch is scarcely modified by the increase of the density stratification, 
which provides an important insight into the global understanding of the impact of 
the density stratification on the stability domain of the dipolar dynamo branch. 
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1 INTRODUCTION 


Observations of low-mass stars reveal very different mag¬ 
netic field topologies, ranging from small-scale fields to large- 
scale dipolar fields, and the last advances in spectropo- 
larimetry should enable one to improve the understanding 


of the magnetic fields of solar-type stars ( Donati &; Land- 
street 2009; Morin et al.|2010 ). Among the three suggestions 
advanced by Larmor to explain the generation of such mag¬ 
netic fields (Larmor 1919), it is now the consensus that their 
decay is prevented by the action of self-excited dynamos in¬ 
duced by the turbulent motions that occur in stellar inte¬ 
riors. More often, these motions are assumed to be driven 
by convection, owing to the temperature difference between 
the inner core and the cooler surface. In dynamo theory, 
this partial transfer of the kinetic energy of a conducting 
fluid into magnetic energy is an instability process: above 
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a certain threshold, electrical currents start to be amplified 
by the fluid flow, so that a magnetic field can be sustained 
against the resistive decay due to ohmic dissipation. 

After Glatzmaier & Roberts (1995), numerical mod¬ 
elling of self-consistent dynamos underwent considerable de¬ 
velopment (in contrast with the small number of success¬ 
ful experimental studies). However, despite the continuous 
increase of computer power, direct numerical simulations 
still face the difficulty to resolve a vast range of spatial 
and temporal scales when attempting to simulate a three- 
dimensional turbulent flow on a magnetic diffusion time- 
scale. As a simplification, one usually resorts to some con¬ 
vective approximations, and most of the early studies were 
relying on the Boussinesq approximation, which performs 
well as long as variations in pressure hardly affect the den¬ 
sity of the fluid. However, this assumption is not valid to 
describe convection in large stratified systems such as stars 
or gas giants, in which the density typically varies over many 
scale-heights between the top and bottom of the convection 
zone. This limitation of the Boussinesq approximation is ba- 
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sically what motivated the use of the anelastic approxima¬ 
tion, originally developed to study atmospheric convection 
( Ogura k Phillips]|1962| Gough||1969 ), to model convection 


in the Earth core and stellar interiors. Indeed, if we assume 
that the overall system remains close to an adiabatically 
stratified reference state at marginal stability so that con¬ 
vective motions can be treated as small perturbations (which 
in turns implies that typical velocities remain small compare 
to the speed of sound), then the anelastic approximation al¬ 
lows us to take some stratification into account while filter¬ 
ing out sound waves for faster numerical integration. This 
approximation can be found in the literature under slightly 


sky k Roberts| 1995;|Lantz Fan 1999; 

Anufriev, Jones k 

Soward 2005; Berkoff, Kersale k Tobias 

2010; Jones et al. 


201 1 | |Alboussiere k Ricard |2013| ), which are in part com¬ 
pared in Brown, Vasil & Zweibel ( 2012|). 


Just as in Boussinesq models (Christensen & Aubert 
2006 |Schrinner, Petitdemange V Dormy|2012| |Yadav et al. 


2013), magnetic fields obtained in anelastic simulations 
^ Gastine, Duarte k Wicht ||2012 Duarte, Gastine k Wicht 


2013 Schrinner et al. 2014) fall into two categories: dipo¬ 


lar dynamos, dominated by a large-scale axial dipole com¬ 
ponent, and multipolar dynamos, characterized by a more 
complex field topology with higher spatial and temporal 
variability. However, these studies identified several differ¬ 
ences specific to anelastic dynamos. For instance, dipolar 
solutions seem more difficult to obtain as the density strat¬ 
ification is increased (Gastine, Duarte k Wicht 2012 Jones 
2014). We found in Schrinner et al. (2014) that for a given 
N e , E and Pr, there seems to exist a critical magnetic 
Pm c below which the dipolar solution is not stable, and 
the higher the density stratification, the higher this criti¬ 
cal magnetic Prandtl number. Furthermore, multipolar dy¬ 
namos with a magnetic field configuration dominated by an 
equatorial dipole seem more easily realized with anelastic 
models than with Boussinesq models. However, we show in 


Raynaud, Petitdemange k Dormy (2014) that this charac¬ 


teristic also stands for weakly stratified models, since it is 
actually related to the use of different mass distributions. In¬ 
deed, the gravity profile may strongly influence the localiza¬ 
tion of the convective cells, depending on whether one con¬ 
siders a homogeneous (g oc r) or a central mass (g oc 1 /r 2 ) 
distribution: as opposed to the former, the latter results in 
the concentration of the convective cells close to the inner 
sphere, which favours the emergence of a less diffusive large- 
scale m— 1 mode at the outer surface of the model. 

Our last study of weakly stratified models with a cen¬ 
tral mass distribution naturally constitutes an appropriate 
reference basis from which a detailed understanding of the 
role of the density stratification in anelastic dynamo mod¬ 
els can be achieved. In this paper, we will primarily focus 
on dipolar dynamos. We aim at clarifying apparent contra¬ 
dictions between previous anelastic studies by investigating 
in more detail the evolution of the stability domain of the 
dipolar branch when increasing the density stratification. To 
that end, we rely on a systematic parameter study consist¬ 
ing of 119 three-dimensional, self-consistent dynamo mod¬ 
els obtained by direct numerical simulations. As opposed to 
previous studies that were focusing on Jupiter’s magnetic 


held ( 

Duarte, Gastine k Wicht 

to 

o 

CO 

Gastine et al. 

2014 

Jones 

2014), we do not consider here 

more realistic mod- 


els to reproduce a particular observation, but instead try 
to understand systematic and general tendencies in anelas¬ 
tic models, as a function of the physical control parameters. 
The anelastic equations are recalled in Section [2] and we 
present our results in Section [3] The complete list of numer¬ 
ical simulations performed for this study is given in Table [AT] 
(see Appendix [Al. 


2 EQUATIONS AND METHODS 


Following Jones et al. (2011), we rely on the LBR formula¬ 
tion of the anelastic approximation ( Lantz k Fan|1999 Bra¬ 


ginsky k Roberts 1995 ). Actually, both the model and the 
numerical methods used here are the same as in ISchrinnerl 


et al. (2014) and Raynaud, Petitdemange k Dormy (2014) 


but we briefly recall them for completeness. We consider a 
spherical shell of width d and aspect ratio x? rotating about 
the z-axis at angular velocity H and filled with a perfect, 
electrically conducting gas with kinematic viscosity is, ther¬ 
mal diffusivity ft, specific heat c p and magnetic diffusivity r] 
(all assumed to be constant). Convection is driven by an im¬ 
posed entropy difference As between the inner and the outer 
boundaries, and the gravity is given by g = -GMr/r 2 , 
where G is the gravitational constant and M the central 
mass. 

The reference state is given by the polytropic equilib¬ 
rium solution of the anelastic system 


P = P c w n+1 , Q=Q c w n , T = T c w, w = co + —, 


Co = 


2wo - X - 1 
1-X ’ 


Cl = 


with 


w 0 = 


X + 1 


(l + x)(l - Wo) 

(i-x ) 2 ’ 


1 + X-Wo 


(1) 

( 2 ) 


(3) 


Xexp(N e /n) + 1’ X 

In the above expressions, n is the polytropic index and 
N g = In {qi/qo) the number of density scale-heights. The 
values P c , g c , and T c are the reference-state density, pres¬ 
sure, and temperature mid-way between the inner and outer 
boundaries, and serve as units for these variables. 

Length is scaled by the shell width d, time by the mag¬ 
netic diffusion time d 2 /r/ and entropy by the imposed en¬ 
tropy difference As. The magnetic field is measured in units 
of y/QQc/JLTj, where /i is the magnetic permeability. Then, the 
equations governing the system are 


d~v 

— + (v • V) v = Pm 


1 P' Pm s ^ 2 ^ 

— — V- 1 - -^Ra — r — — z x v 


E 


+ + 


E w 1 


Pr 

-(V x B) x B 


E 


f = Vx(vxB) + V»B, 


ds —n—i Pm 

— +vVs = w - 5 -V 

at Pr 


—V • (w n+ Vs) 


+ — [E- 
w L 


1 —n 
W 


(V x B) 2 + Q v 


V • (w n v) = 0 , 
V • B = 0 . 


(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 
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Figure 1 . Left: dipolar (black circles) and multipolar (white squares) dynamos as a function of Ra/Ra c and Pm, for N e = 0.5 (a), 
N g = 1.5 (c) and N g = 2.0 (e). A cross indicates the absence of a self-sustained dynamo. Right: the relative axial dipole field strength 
/dipax versus the local Rossby number for N g = 0.5 (b), N g = 1.5 (d) and N g = 2.0 (f). 
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The viscous force F*, in Eq. @ is given by F v — w n VS, 
where S is the rate of strain tensor 


f u ^ \ 

1 . 

( dvi 

d y j \ 

( - ^OijV • V 1 , 

Gij ~~ 2 ^ 


+ d^J 


(9) 

Moreover, the expressions of the dissipation parameter Di 
and the viscous heating Q u in Eq. © are 

c\Pr 


Di = 


PmRa 


and 


Qu = 2 


etjeij - -(V • v) s 


( 10 ) 

( 11 ) 


We impose stress-free boundary conditions for the ve¬ 
locity field at both the inner and the outer spheres, the mag¬ 
netic held matches a potential held inside and outside the 
huid shell, and the entropy is hxed at the inner and outer 
boundaries. Besides, both weak and strong held initial con¬ 
ditions have been tested for all models, since the system 
may exhibit hysteretic transitions between dynamo branches 
when stress-free boundary conditions are used. 

The system of equations ©-© involves seven con¬ 
trol parameters, namely the Rayleigh number Ra = 
GMdks/ (vkc v ), the Ekman number E — z^/(Qd 2 ), the 
Prandtl number Pr = v / k, and the magnetic Prandtl num¬ 
ber Pm = vI 77 , together with the aspect ratio x, the poly¬ 
tropic index n, and the number of density scale-heights N g 
that dehne the reference state. We choose E = 10 —4 , Pr = 1 


and n — 2 . Different from Gastine, Duarte & Wicht ( 2012 ), 


we also kept the central gravity prohle and the aspect ratio 
X = 0.35 hxed for all simulations, but varied the magnetic 
Prandtl number, which turns out to be a key point to un¬ 
derstand the partial divergence of our conclusions. 

The equations are integrated in average for one mag- 
netic diffusion time with the anelastic version of p arody 
( Dormy, Cardin &; Jault||l998 Schrinner et al.|2014| )p|The 

vector helds are transformed into scalars using the poloidal- 
toroidal decomposition. The equations are then discretized 
in the radial direction with a finite-difference scheme; on 
each concentric sphere, variables are expanded using a spher¬ 
ical harmonic basis. The coefficients of the expansion are 
identihed with their degree £ and order m. Typical resolu¬ 
tions are 288 points in the radial direction (up to 320 points). 
The spectral decomposition is truncated at a hundred modes 
(up to <£ max ~ ra max < 128), in order to observe for both 
spectra a decrease of more than two orders of magnitude 
over the range of l and m. The highest resolutions are re¬ 
quired for the models with the highest density stratification 

(•V, - 3). 

The amplitudes of the velocity and the magnetic fields 
are measured in terms of the Rossby number Ro = 
\f2E^E/Pm and Lorentz number Lo = \/2E m E / Pm, where 
Ek and E m are the energy densities integrated over the fluid 
shell, 


E k = - 

2E 


L 


w n v 2 dv and 


P 1 Pm f 2 . 

e ‘-2v-etJ v b *’■ 

( 12 ) 


1 The integration times range from 0.63 to 5.2 magnetic diffusion 
times (for the models 101m and 004m, respectively). 


Likewise, the measure of the mean zonal flow is given by the 
zonal Rossby number Ro z based on the averaged toroidal 
axisymmetric kinetic energy. 

We also define a local Rossby number Rot — Ro c £ c /ty 
based on the mean harmonic degree £ c of the velocity com¬ 
ponent v c from which the mean zonal flow has been sub¬ 
tracted (Schrinner, Petitdeman ge &; Dormy|2012 Schrinner 


et al.[|2Q14D, 


= £ 


, (w n (vc)r(vc)t) 
(w n V c • V c ) 


(13) 


where the brackets denote an average over time and radii. 
The contribution of the mean zonal flow is removed for cal¬ 
culating Ro c . 

Furthermore, as the stratification is increased, it turns 
out that it is useful to examine the variations with depth 
of the local Rossby number, defined in such a way that 
Rot — f^° Rot(r)r 2 dr. We found that it is more suitable 
to slightly adapt our initial definition and investigate the 
radial dependence of Ro\ (r), which differs from Rot (r) in so 
far as the velocity is not weighted by the reference density 
profile w n . We also checked that, in our range of N g , both es¬ 
timates of a characteristic velocity do not make a qualitative 
difference on the volume-averaged quantities. For instance, 
the difference between the values of the magnetic Reynolds 
number Rm = RoPm/E is about 1 per cent at N e = 0.5. 
Of course, it increases with N e : energy-based estimates lead 
to lower values about 7 and 10 per cent for N g = 2.0 and 
N g = 2.5, respectively. However, this does not change our 
conclusions, and that is why we do not adapt our definition 
for volume-averaged quantities. 


3 RESULTS 
3.1 Bistability 


Figure ^ shows the distribution of dipolar and multipo¬ 
lar dynamos in the parameter space ( Ra/Ra c , Pm) (left- 
hand panels), together with the corresponding dipolarities 
(right-hand panels), for increasing density stratifications 
from top to bottom. One can see that several examples of 
bistable pairs are displayed. Bistability is commonly known 
for Boussinesq and anelastic models, and is related to the use 
of stress-free boundary conditions that allows for the growth 
of stronger zonal winds ( Sasaki et al.||2011| |Schrinner, Pe¬ 


titdemange &; Dormy |2012 Gastine, Duarte &; Wicht| 2012 ) 
For N g = 0.5, the regime diagram in Fig. |l(a)| does not qual¬ 
itatively differ from what we can observe in the Boussinesq 
regime. As we found in |Raynaud, Petitdemange fc Dormy | 
(2014), the multipolar branch undergoes a supercritical bi¬ 


furcation as Ra is increased, whereas the dipolar one still 
loses its stability in favour of the multipolar branch at low 
Rayleigh and magnetic Prandtl numbers. When increasing 
the density contrast N g to 1.5, one can note in Fig. |l(c) that 
the overlap between the two branches shrinks. At N g = 2 , 
we do not observe a bistable case. More generally, for all 
models of our sample with a density stratification N g > 2 , 
the saturated field of the dynamo is not anymore sensitive 
to the amplitude of the initial magnetic field. To understand 
this evolution from Boussinesq models to anelastic models 
with moderate stratification ( N g < 1.5), it is worth stress¬ 
ing that the transition process from the multipolar to the 
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Figure 2. (a): Evolution of the zonal Rossby number as a function of Pm for a dynamo models with Ra = 4x 10 6 , at N g = 0.5 (blue) 
and N g = 1.5 (red). Circles (squares) stands for dipolar (multipolar) dynamos, (b): Dipolar dynamos in the parameter space ( Ra/Ra c , 
Pm), for increasing density stratifications: N g = 0.1 (grey), N g = 0.5 (blue) N g = 1.5 (red), N g = 2.0 (green) and N g = 2.5 (black). 



(a) 



(b) 


Figure 3. (a): Magnetic Reynolds number as a function of N e for dipolar dynamos, (b): Our sample of dipolar (circles) and multipolar 
(squares) dynamos in the parameter space ( Ra/Ra c , N g ). 


dipolar branch triggered by the increase of Pm still applies 


to our sample of models (see Schrinner, Petitdemange & 
D ormy| 2012). Figure [2(a)] illustrates on a few cases the pro¬ 


gressive merging of the multipolar branch which is indeed 
lost when its zonal Rossby number becomes comparable to 
the zonal Rossby number of the dipolar branch. For a given 
Rayleigh number, the fact that the mean zonal flow of the 
multipolar branch decreases with Pm (and eventually be¬ 
comes too small to prevent the growth of the axial dipole) is 
actually the limiting factor of the upper extent of the mul¬ 
tipolar branch in the left-hand panels of Fig. [l] This also 
emphasizes the essential role played by differential rotation 
in the dynamo mechanism of the multipolar branch, often 
accounted for in terms of Q-effect. 

Interestingly, the zonal Rossby number for multipolar 
dynamos substantially decreases between N g = 0.5 and 1.5 


(see the blue and red squares in Fig. |2(a) ), while it remains 
of the same order for dipolar dynamos. Hence, the available 
range of Pm for the multipolar solution is reduced, which 
therefore explains the relative shrinking of the bistable re¬ 
gion when comparing Figs |l(a)| and 1(c) The simplest ar¬ 
gument to understand this downtrend is given by the com¬ 
parison of the x-axis in Fig. ^ which reveals that the dy¬ 
namo onset moves closer to the onset of convection when the 
density stratification is increased, as mentioned by |Gastine,| 


Duarte & Wicht (2012). Indeed, despite changing the value 


of Ng , we found that the Rayleigh numbers we had to con¬ 
sider always stay of the order of 10 6 . At the same time, the 
critical Rayleigh number for the linear onset of convection 
monotonically increases with N e . From table B.l in |Schrim| 


ner et al. (2014), we have in our case the following values of 
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3.34 x 10 5 , 9.25 x 10 5 and finally 1.43 x 10 6 for the sequence 
of density stratifications N e = 0.5, 1.5 and 2.0, respectively. 


3.2 Dipole onset 

The density stratification strongly impacts on the stability 


domain of the dipolar branch, as we clearly see in Fig. 2(b) 


In this figure, we included data from Raynaud, Petitde- 
mange & Dormy (2014) in order to better highlight the dif¬ 


ferences with Boussinesq simulations. For moderate values 
of N g at a fixed Pm , the critical value of Ra/Ra c at which it 
is possible to sustain a dipolar dynamo rapidly falls off (up 
to a factor of 4 if we consider the line Pm = 1). However, 
this tendency hardly persists once we reach Ra/Ra c ~ 5 
for N g = 1.5, and the further increase of N e mainly affects 
the critical magnetic Prandtl number Pm c below which it is 
not possible to sustain a dipolar dynamo. In our sample of 
models, the increase of Pm c becomes effective for N e > 2, 
but we already reported it as a general tendency in |Schrin-| 
ner et al. (2014). Figure [1(c) | enables us to conclude that 
0.5 < Pm c < 0.75 for N e = 1.5, whereas from Fig. 1(e) it 
is clear that Pm c > 1 for N g — 2. 

The fact that dipolar dynamos are found closer to the 
convection threshold as N g increases can be more or less 
readily understood if one notices that, despite the increase 
of the density stratification, the critical magnetic Reynolds 
number Rm c of the dipolar branch does not significantly 
vary, but stays in first approximation of the order of 10 2 , as 
shown in Fig. 3(a)| Then, if we take this as a necessary con¬ 
dition to obtain a dipolar solution, and given the fact that 
for a constant value of Ra/Ra c the flow amplitude increases 
with N g (Gastine, D uarte &; Wicht||2012 ), it explains why 
the dipolar branch can be found closer to the onset of con¬ 
vection when the stratification increases. However, we will 
see in the next subsection that, as N g is further increased, 
not only does the dipolar branch occur closer to the onset of 
convection, but also higher magnetic Prandtl numbers have 
to be considered to maintain a sufficiently high Rm while 
preventing the collapse of the dipole. 


3.3 Dipole collapse 


Another striking feature that arises when investigating the 
stability domain of the dipolar branch is that the range of 
Rayleigh numbers over which it extends becomes smaller 
and smaller as N g is increased. This is clearly visible in 
Fig. |3(b) that shows for different N g the transition from 
the dipolar to the multipolar branch resulting from the in¬ 
crease of Ra. In other words, at this moderate value of the 
Ekman number, dipolar dynamos are confined in a narrower 
and narrower window of Rayleigh numbers, which explains 
why dipolar solutions may seem more difficult to obtain at 
higher N e , despite comparable critical magnetic Reynolds 
numbers. As for the modification of the dynamo onset, this 
can be related to the fact that for a given value of Ra/Ra c , 
the Rossby number Ro increases with A^^The transition 



Figure 4. Average values of /dip ax f° r dipolar dynamos as a 
function of N g . Error bars represent the standard deviation. The 
average is done with 11 models for N g = 2.5. 


from a dipolar to a multipolar solution triggered by an in¬ 
crease of Ra is related to the fact that inertia becomes sig¬ 
nificant in the force balance. We know from IChristensen fel 


Aubert (2006) that this transition can be measured by a 


local Rossby number Roi based on a characteristic length- 
scale of the flow. We find that the collapse of the dipole 
still occurs for Roi 


0.1 when N g < 2, which is consistent 


with the results in Gastine, Duarte &; Wicht (2012). This 
is illustrated by Figures [l(b)| 1 (d) | and |l(f)| which show the 
relative axial dipole field strength /di Pax computed at the 
outer sphere, as a function of Ro£. In Fig. l(b)| the very low 
values of /dip ax at low Roi are characteristics of multipo¬ 
lar dynamos dominated by an equatorial dipole component. 
We showed in Raynaud, Pet it demange & Dormy (2014) that 
this magnetic configuration arises close to the dynamo on¬ 
set and when convective cells are localized close to the in¬ 
ner sphere. However, we know from hydrodynamic studies 
that the convection cells move towards the outer shell when 
the stratification is increased ( |Jones, Kuzanyan &; Mit chell 


2009; Gastine fc Wicht||2012 ), which explains why this fea¬ 


ture tends to disappear in Figs |l(d~ and 1(f) Besides, we 
see in Fig. [4] that the values of /di Pax tend to decrease with 
N e , which is also clear if we focus for instance on the dipolar 
branch in Fig. |1 (f )] for which /di Pax < 0.8. As expected, this 
indicates that the small magnetic scales at the outer surface 
are favoured with the increase of the stratification. This is 
also clearly confirmed by the comparison of the radial mag¬ 
netic fields at the outer surface of the model, as shown in 
the left-hand panelsleft of Figs [5] and [6] Finally, we also re¬ 
port the existence of multipolar dynamos whose dipolarity 
displays strong variations in time. This leads to averaged 
values of /di Pax ~ 0.5, as one can notice in Fig. 1(d) These 
dynamos usually exhibit a relatively strong axial dipole com¬ 
ponent which undergoes reversals during which the value of 
fd i Pax decreases drastically. Duarte, Gastine & Wicht (2013) 
also reported similar behaviour for dynamo models with a 
variable electrical conductivity. 


2 For instance, one can compare the models 008d, 055d and 083d 
for which Ra/Ra c ~ 6 and an increasing Ro of 4.8 X 10 3 , 1.3 X 
10 -2 and 1.8 X 10 2 , respectively, or else the models 021d, 051d 


and 096d that have a similar Rossby number of 0.01, but for which 
Ra/Ra c is about 12, 5 and 2.9, respectively. 
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(a) N e = 2.5 (b) N q = 2.5 

Figure 6. Snapshot of B r (r = r Q ) (a) and equatorial cut of v r (b) for a dipolar dynamo with N g = 2.5, Pm = 2, Ra = 7.40xl0 6 = 3ARa c . 


For N g > 2.0, we found that the dipole collapse tends 
to occur at values of Ron lower than 0.1. However, it is likely 
that a volume-averaged quantity becomes less relevant when 
applied to models with a substantial stratification. For in¬ 
stance, we see in Figs |5(b)| and [6(b)] that the smaller struc¬ 
tures that develop at N e = 2.5 are confined close to the 
outer boundary, whereas there are no significant differences 
in the radial flow at mid-depth. Thus, we also examined the 
radial dependence of the different components of the local 
Rossby number Ro*, which is computed as the product of 
two terms: a convective Rossby number based on the velocity 
field v c from which the mean zonal flow has been subtracted 
(see Fig. 7(a)) and a characteristic length-scale based on the 
mean harmonic degree of v c (see Fig. |7(b)| ). We find that 
the monotonicity of Ro* changes as N g is increased. Indeed, 


for low stratifications, Ro*(r) mainly decreases with radius, 
whereas for N g > 2.5 it becomes an increasing function of r 
that steepens slightly close to the outer surface. Figure |7(cj| 
shows the evolution of Ro*(r) for increasing Rayleigh num¬ 
bers up to the loss of the dipolar solution, at N g = 2.5 and 
3.0. When the transition to the multipolar branch is reached, 
we see that Rol tends to increase faster close to the outer 
surface, while the volume-averaged value can stay below the 
critical value of 0.1. Thus, it seems that inertia still causes 
the collapse of the dipolar branch, despite the fact that the 
usual local Rossby number criterion is not appropriate to 
separate the two dynamo branches for significant density 
stratifications. 
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Figure 7. The convective Rossby number (a), the convective length-scale (b) and the local Rossby number (c) as a function of radius 
for dipolar (solid lines) and multipolar (dashed lines) dynamos at (N g = 2.5, Pm = 2) (thin lines) and (N g = 3, Pm = 4) (thick lines). 



3.4 Dynamo mechanisms 

Finally, we try to investigate whether the dynamo mech¬ 
anisms at work on the dipolar branch are modified when 
the stratification is increased. We see in Fig. [8] that the axi- 
symmetric azimuthal magnetic field we observe at N g — 1.5 
is strongly reminiscent of the magnetic structures that can 
be observed with Boussinesq models, which are u sually in- 
terpreted in terms of o? dynamos (Olson, Christensen h 
Glatzmaier||1999] ISchrinner et al7| |2007| |Schrinner, Petit- 
demange & Dormy||2011| |2012| ). tnside the tangent cylin¬ 
der, the azimuthal magnetic field is mainly produced by the 
Q-effect, rB r <9 (r _1 \©) /<9r+r _1 sinO'Bed (sin0 -1 V^) / dO , 
which correlates inside the tangent cylinder with the axisym- 
metric azimuthal magnetic field, when comparing Figs |8(b~ 
and 8 (c) | However, outside the tangent cylinder, the most 
part of the mean azimuthal field does not seem to be the re¬ 
sult of the 12-effect, and it is thus likely that the essential re¬ 
generation of the poloidal field is achieved by a-effect, lead¬ 
ing to the emergence of characteristic equatorial patches of 
opposite polarity (see e.g. Christensen 2011 Schrinner, Pe- 


titdemange & Dorm y|2012 ). We did not find in our sample of 
models tangible evidence that would invalidate this scenario 
at higher N g . For instance, at N g = 2.5, we see in Fig. [9] 
that the major differences lie in the stronger axisymmetric 
azimuthal velocity (compare Figs |8(a)] and |9(a)] ). Neverthe¬ 
less, the axisymmetric azimuthal magnetic field shown in 
Fig. |9(b) seems only modified about a colatitude 0 ^ yr/4 
close to the outer surface, and keeps now the same polarity in 
each hemisphere outside the tangent cylinder. This change 
can be correlated to the modifications of the axisymmetric 
azimuthal velocity, which in turn affect the 12-effect (com¬ 
pare Figs |8(cJ and 9(c) | ) . Of course, we are for now limited 
to the observation of correlations, but it would be interest¬ 
ing to have a further insight into the dynamo mechanism 
in anelastic simulations using a test field method, in the 
spirit of the Boussinesq study by |Schrinner, Petitdemange| 


& Dormy (2012). 


© 2014 RAS, MNRAS 000, [1^13] 
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4 CONCLUSION 

With this systematic parameter study, we are able to im¬ 
prove our understanding of the successive modifications that 
are exhibited by the stability domain of the dipolar branch 
when increasing the density stratification in anelastic dy¬ 
namo models. In general, dipolar dynamos are found closer 
to the onset of convection. Moreover, we show that dipole- 
dominated solutions can be observed even at high density 
stratifications, provided high enough magnetic Prandtl num¬ 
bers are considered. Besides, this study also highlights why 
dipolar dynamos seem more difficult to find in anelastic sim¬ 
ulations. Indeed, this tendency is usually reported as a gen¬ 
eral statement, but here we show that this impression mainly 
results from the fact that the dipolar branch extends on a 
smaller and smaller range of Rayleigh numbers as N e is in¬ 
creased. However, despite the relative shrinking of the sta¬ 
bility domain, we found that the critical magnetic Reynolds 
number of the dipolar branch seems scarcely modified in 
the overall process. At the same time, the higher N e , the 
faster convection will develop as we depart from the on¬ 
set. In consequence, the higher N e , the faster is reached 
the critical Rossby number above which inertia causes the 
collapse of the dipole. This explains why dipolar dynamos 
become clearly confined in a smaller region of the param¬ 
eter space. However, we stress that, in terms of magnetic 
Reynolds number, the dynamo threshold does not signifi¬ 
cantly increase with the density stratification in the range 
of N g we investigated. 

In addition, this study also suggests that the scarcity of 
dipolar solutions for substantial density stratifications would 
thus rather come from the restriction of the parameter space 
being currently explored (because of computational limita¬ 
tions), rather than an intrinsic modification of the dynamo 
mechanisms that would be caused by the density stratifica¬ 
tion. Furthermore, if we decrease the Ekman number from 
E = 10 -4 to 3 x 10 -5 keeping Pr = 1, we find that we re¬ 
cover three examples of bistable pairs at N e = 2, for Pm = 1 


at Ra/Ra c = 2.6 and for Pm G {1,2} at Ra/Ra c = 2.9. 
Then, beyond the results of this study, and for low values 
of the Ekman number that are currently very expensive to 
simulate, it seems more likely that dipolar solutions will per¬ 
sist in a larger region of parameter space (see also 
20T4||Jones|2 01~). 


Duarte 


Despite the fact that it is not straightforward to re¬ 


late the output of numerical models with observations (Gas- 
tine et al.|2013 ), the bistability that is reported for numeri¬ 


cal simulations can be similarly observed with real objects. 
For instance, in a spectropolarimetric survey done with a 
sample of active M dwarfs, Morin et al. (2010) report two 
distinct categories of magnetic topologies. They distinguish 
strong axisymmetric dipolar fields and weak fields with sig¬ 
nificant non-axisymmetric components, and both configura¬ 
tions seem to be observed on objects with similar stellar pa¬ 
rameters. After Schrinner, Petitdemange & Dormy (2012), 
we show that the bistable behaviour observed in numerical 
models could be a possible way towards a better understand¬ 
ing of the broad diversity of the magnetic fields of M dwarfs, 
and that it cannot be ruled out even when taking into ac¬ 
count the density stratification. The understanding of the 
impact of the stratification on the dynamo mechanisms de¬ 
serves further studies. 
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Table Al: Overview of the simulations carried out, with E — 10 4 
Pr = 1, x — 0.35, and n — 2. 


model 

N e 

Ra 


Pm 


Ro 


Ron 


Ro z 


Lo 

/ di Pax 

001m 

0.5 

1.500 

X 

10 s 

2.00 

3.2 

X 

nr® 

2.0 

X 

IO -2 

1.8 

X 

nr 3 

2.3 

X 

10“ 3 

6.6 

X 

10“ 2 

002m 

0.5 

1.750 

X 

10 6 

1.00 

3.7 

X 

10“ 3 

2.2 

X 

10“ 2 

2.7 

X 

10“ 3 

1.8 

X 

10“ 3 

1.6 

X 

10 -2 

003m 

0.5 

1.800 

X 

10® 

0.75 

3.5 

X 

10“ 3 

2.8 

X 

10“ 2 

2.8 

X 

IO -3 

1.9 

X 

10“ 3 

2.2 

X 

10 -3 

004m 

0.5 

1.850 

X 

10® 

1.00 

3.9 

X 

nr 3 

2.4 

X 

IO -2 

4.2 

X 

nr 3 

2.0 

X 

10“ 3 

2.5 

X 

10“ 2 

005m 

0.5 

2.000 

X 

10® 

0.75 

4.1 

X 

10“ 3 

2.3 

X 

10“ 2 

5.2 

X 

10“ 3 

2.4 

X 

10“ 3 

9.0 

X 

10- 5 

006m 

0.5 

2.000 

X 

10® 

1.50 

5.5 

X 

10“ 3 

2.8 

X 

10“ 2 

3.0 

X 

IO -3 

3.1 

X 

10“ 3 

5.3 

X 

10“ 2 

007m 

0.5 

2.000 

X 

10® 

2.00 

5.3 

X 

nr 3 

2.7 

X 

10“ 2 

2.0 

X 

nr 3 

3.7 

X 

10“ 3 

9.4 

X 

10“ 2 

008m 

0.5 

2.000 

X 

10® 

3.00 

5.1 

X 

10“ 3 

2.7 

X 

10“ 2 

1.7 

X 

10“ 3 

4.1 

X 

10“ 3 

2.7 

X 

10- 1 

008d 

0.5 

2.000 

X 

10® 

3.00 

4.8 

X 

10“ 3 

2.6 

X 

10“ 2 

1.0 

X 

IO -3 

5.5 

X 

10“ 3 

7.7 

X 

10 _1 

009d 

0.5 

2.000 

X 

10® 

5.00 

4.8 

X 

nr 3 

2.5 

X 

10“ 2 

1.1 

X 

nr 3 

5.7 

X 

10“ 3 

6.8 

X 

10- 1 

010m 

0.5 

2.500 

X 

10® 

0.75 

7.7 

X 

10“ 3 

3.3 

X 

10“ 2 

5.2 

X 

10“ 3 

3.8 

X 

10“ 3 

5.1 

X 

10 -2 

011m 

0.5 

2.500 

X 

10® 

1.00 

7.2 

X 

10“ 3 

3.5 

X 

10“ 2 

3.6 

X 

IO -3 

4.5 

X 

10“ 3 

7.6 

X 

10“ 2 

012m 

0.5 

2.500 

X 

10® 

1.50 

7.1 

X 

nr 3 

3.5 

X 

10“ 2 

3.1 

X 

nr 3 

4.5 

X 

10“ 3 

3.3 

X 

io- 1 

013d 

0.5 

2.500 

X 

10® 

4.00 

6.0 

X 

10“ 3 

3.3 

X 

10“ 2 

1.5 

X 

10“ 3 

6.9 

X 

10“ 3 

7.3 

X 

lO -1 

014m 

0.5 

3.000 

X 

10® 

1.00 

8.6 

X 

10“ 3 

4.2 

X 

10“ 2 

4.0 

X 

IO -3 

5.5 

X 

10“ 3 

1.3 

X 

10 _1 

015m 

0.5 

3.000 

X 

10® 

2.00 

8.2 

X 

nr 3 

4.2 

X 

10“ 2 

2.8 

X 

nr 3 

6.3 

X 

10“ 3 

3.0 

X 

io - 1 

015d 

0.5 

3.000 

X 

10® 

2.00 

7.5 

X 

10“ 3 

3.9 

X 

10“ 2 

1.7 

X 

10“ 3 

7.9 

X 

10“ 3 

7.5 

X 

IO -1 

016d 

0.5 

3.000 

X 

10® 

3.00 

7.5 

X 

10“ 3 

4.0 

X 

10“ 2 

1.7 

X 

IO -3 

8.5 

X 

10“ 3 

7.2 

X 

IO -1 

017d 

0.5 

3.000 

X 

10® 

4.00 

7.5 

X 

nr 3 

4.0 

X 

10“ 2 

1.7 

X 

nr 3 

8.5 

X 

10“ 3 

7.2 

X 

IO- 1 

018d 

0.5 

3.000 

X 

10® 

5.00 

7.5 

X 

10“ 3 

4.1 

X 

10“ 2 

1.8 

X 

10“ 3 

8.8 

X 

10“ 3 

6.8 

X 

IO- 1 

019d 

0.5 

3.000 

X 

10® 

6.00 

7.6 

X 

10“ 3 

4.1 

X 

10“ 2 

1.9 

X 

IO -3 

9.0 

X 

10“ 3 

6.4 

X 

IO -1 

020m 

0.5 

4.000 

X 

10® 

0.50 

1.3 

X 

10“ 2 

5.3 

X 

10“ 2 

9.3 

X 

10“ 3 

6.6 

X 

10“ 3 

2.3 

X 

IO- 1 

021m 

0.5 

4.000 

X 

10® 

1.00 

1.2 

X 

10“ 2 

5.6 

X 

10“ 2 

5.7 

X 

10“ 3 

7.9 

X 

10“ 3 

2.8 

X 

IO- 1 

021d 

0.5 

4.000 

X 

10® 

1.00 

1.0 

X 

10“ 2 

5.4 

X 

10“ 2 

2.0 

X 

10 -3 

1.1 

X 

10“ 2 

8.7 

X 

IO -1 

022m 

0.5 

4.000 

X 

10® 

2.00 

1.5 

X 

10“ 2 

5.7 

X 

10“ 2 

3.5 

X 

10“ 3 

8.8 

X 

10“ 3 

2.6 

X 

IO- 1 

022d 

0.5 

4.000 

X 

10® 

2.00 

1.0 

X 

10“ 2 

5.4 

X 

10“ 2 

1.9 

X 

10“ 3 

1.1 

X 

10“ 2 

6.6 

X 

IO- 1 

023d 

0.5 

4.000 

X 

10® 

3.00 

1.0 

X 

10“ 2 

5.3 

X 

10“ 2 

2.0 

X 

10 -3 

1.1 

X 

10“ 2 

6.8 

X 

10 _1 

024d 

0.5 

4.000 

X 

10® 

4.00 

1.0 

X 

10“ 2 

5.3 

X 

10“ 2 

1.9 

X 

10“ 3 

1.1 

X 

10“ 2 

6.5 

X 

IO- 1 

025d 

0.5 

4.000 

X 

10® 

6.00 

1.0 

X 

10“ 2 

5.3 

X 

10“ 2 

2.2 

X 

10“ 3 

1.2 

X 

10“ 2 

5.4 

X 

IO- 1 

026m 

0.5 

5.000 

X 

10® 

0.50 

1.5 

X 

10“ 2 

6.8 

X 

10“ 2 

8.4 

X 

10 -3 

8.8 

X 

10“ 3 

2.3 

X 

10 _1 

027m 

0.5 

5.000 

X 

10® 

1.00 

1.4 

X 

10“ 2 

7.1 

X 

10“ 2 

5.5 

X 

10“ 3 

1.0 

X 

10“ 2 

2.5 

X 

IO- 1 

027d 

0.5 

5.000 

X 

10® 

1.00 

1.3 

X 

10 -2 

6.7 

X 

10“ 2 

2.0 

X 

10 -3 

1.4 

X 

10“ 2 

8.1 

X 

10 _1 

028d 

0.5 

5.000 

X 

10® 

2.00 

1.3 

X 

10“ 2 

6.6 

X 

10“ 2 

2.2 

X 

10 -3 

1.4 

X 

10“ 2 

6.8 

X 

10 _1 

029d 

0.5 

5.000 

X 

10® 

3.00 

1.3 

X 

10“ 2 

6.8 

X 

10“ 2 

2.1 

X 

10“ 3 

1.4 

X 

10“ 2 

7.2 

X 

IO- 1 

030d 

0.5 

5.000 

X 

10® 

4.00 

1.3 

X 

IO -2 

6.5 

X 

10“ 2 

2.5 

X 

10 -3 

1.6 

X 

10“ 2 

7.0 

X 

10 _1 

031d 

0.5 

5.000 

X 

10® 

5.00 

1.3 

X 

10“ 2 

6.7 

X 

10“ 2 

2.5 

X 

10 -3 

1.5 

X 

10“ 2 

6.8 

X 

IO -1 

032m 

0.5 

6.000 

X 

10® 

1.00 

1.7 

X 

10“ 2 

8.5 

X 

10“ 2 

6.3 

X 

10“ 3 

1.3 

X 

10“ 2 

2.8 

X 

IO- 1 

032d 

0.5 

6.000 

X 

10® 

1.00 

1.5 

X 

IO -2 

8.0 

X 

10“ 2 

2.6 

X 

10 -3 

1.7 

X 

10“ 2 

8.2 

X 

10 _1 

033m 

0.5 

7.000 

X 

10® 

1.00 

2.0 

X 

10“ 2 

9.6 

X 

10“ 2 

7.4 

X 

10 -3 

1.4 

X 

10“ 2 

2.5 

X 

10 _1 

033d 

0.5 

7.000 

X 

10® 

1.00 

1.8 

X 

10“ 2 

8.6 

X 

10“ 2 

3.0 

X 

10“ 3 

2.0 

X 

10“ 2 

8.5 

X 

IO- 1 

034m 

0.5 

9.000 

X 

10® 

1.00 

2.5 

X 

IO -2 

1.2 

X 

IO -1 

8.2 

X 

10 -3 

1.9 

X 

10“ 2 

3.6 

X 

IO -1 

035m 

0.5 

1.000 

X 

10 7 

1.00 

2.7 

X 

10“ 2 

1.3 

X 

IO -1 

8.7 

X 

10 -3 

2.1 

X 

10“ 2 

3.4 

X 

10 _1 

036m 

1.5 

2.500 

X 

10® 

0.75 

3.9 

X 

10“ 3 

2.6 

X 

10“ 2 

2.9 

X 

10“ 3 

3.3 

X 

10“ 3 

4.7 

X 

IO -2 

037m 

1.5 

2.500 

X 

10® 

1.00 

4.5 

X 

IO -3 

2.7 

X 

10“ 2 

2.0 

X 

10 -3 

3.6 

X 

10 -3 

1.9 

X 

10 _1 

038m 

1.5 

2.500 

X 

10® 

1.50 

4.3 

X 

IO -3 

2.8 

X 

10“ 2 

9.2 

X 

10 -4 

4.4 

X 

10“ 3 

1.0 

X 

10 _1 

038d 

1.5 

2.500 

X 

10® 

1.50 

4.3 

X 

10“ 3 

2.8 

X 

10“ 2 

8.2 

X 

10“ 4 

5.0 

X 

10“ 3 

7.9 

X 

IO- 1 

039d 

1.5 

2.500 

X 

10® 

2.00 

4.3 

X 

IO -3 

3.0 

X 

10“ 2 

7.7 

X 

10 -4 

4.6 

X 

10 -3 

7.1 

X 

10 _1 

040d 

1.5 

2.500 

X 

10® 

3.00 

4.2 

X 

IO -3 

2.9 

X 

10“ 2 

8.0 

X 

10 -4 

4.8 

X 

10“ 3 

7.3 

X 

10 _1 

041m 

1.5 

3.000 

X 

10® 

0.75 

5.5 

X 

10“ 3 

3.6 

X 

10“ 2 

2.9 

X 

10“ 3 

4.8 

X 

10“ 3 

2.1 

X 

IO- 1 

042m 

1.5 

3.000 

X 

10® 

1.00 

6.0 

X 

IO -3 

3.6 

X 

10“ 2 

2.0 

X 

10 -3 

5.1 

X 

10 -3 

3.9 

X 

10 _1 

042d 

1.5 

3.000 

X 

10® 

1.00 

5.7 

X 

IO -3 

3.7 

X 

10“ 2 

1.1 

X 

10 -3 

6.4 

X 

10“ 3 

8.5 

X 

IO -1 

043d 

1.5 

3.000 

X 

10® 

2.00 

5.5 

X 

10“ 3 

3.7 

X 

10“ 2 

1.1 

X 

10“ 3 

6.2 

X 

10“ 3 

7.1 

X 

IO- 1 

044d 

1.5 

3.700 

X 

10® 

3.00 

7.3 

X 

IO -3 

4.7 

X 

10“ 2 

1.4 

X 

10 -3 

7.9 

X 

10 -3 

6.0 

X 

10 _1 

045m 

1.5 

4.000 

X 

10® 

0.50 

9.5 

X 

IO -3 

4.9 

X 

10“ 2 

5.3 

X 

10 -3 

6.2 

X 

10“ 3 

2.8 

X 

10 _1 

046m 

1.5 

4.000 

X 

10® 

0.75 

8.7 

X 

10“ 3 

5.2 

X 

10“ 2 

3.0 

X 

10“ 3 

6.9 

X 

10“ 3 

4.5 

X 

IO- 1 

046d 

1.5 

4.000 

X 

10® 

0.75 

8.5 

X 

IO -3 

5.3 

X 

10“ 2 

1.5 

X 

10 -3 

9.3 

X 

10 -3 

8.2 

X 

IO -1 

047d 

1.5 

4.000 

X 

10® 

1.00 

8.5 

X 

IO -3 

5.2 

X 

10“ 2 

1.4 

X 

10 -3 

9.4 

X 

10“ 3 

7.8 

X 

10 1 
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12 


Table Al: (continued) 


model 

N e 

Ra 


Pm 


Ro 


Roe 


Ro z 


Lo 

/diPax 

048d 

1.5 

4.000 

X 

10® 

2.00 

8.5 

X 

10 -3 

5.3 

X 

10“ 2 

1.5 

X 

10 -3 

8.4 

X 

KT 3 

6.9 

X 

lO -1 

049m 

1.5 

4.625 

X 

10® 

0.50 

9.5 

X 

10 -3 

5.8 

X 

10“ 2 

6.0 

X 

10 -3 

7.5 

X 

10“ 3 

4.5 

X 

10 _1 

050m 

1.5 

4.625 

X 

10 6 

0.75 

1.0 

X 

10“ 2 

6.2 

X 

10“ 2 

3.6 

X 

10“ 3 

8.1 

X 

10“ 3 

5.3 

X 

10- 1 

050d 

1.5 

4.625 

X 

10® 

0.75 

1.0 

X 

10“ 2 

6.2 

X 

10“ 2 

1.7 

X 

10 -3 

1.1 

X 

10“ 2 

8.6 

X 

10 _1 

051d 

1.5 

4.625 

X 

10® 

1.00 

1.0 

X 

10“ 2 

6.3 

X 

10“ 2 

1.7 

X 

10 -3 

1.1 

X 

10“ 2 

7.6 

X 

10 _1 

052d 

1.5 

5.000 

X 

10® 

1.00 

1.1 

X 

10“ 2 

6.7 

X 

10“ 2 

1.8 

X 

10“ 3 

1.2 

X 

10“ 2 

7.5 

X 

10- 1 

053d 

1.5 

5.000 

X 

10® 

2.00 

1.1 

X 

10“ 2 

6.7 

X 

10“ 2 

1.8 

X 

10 -3 

1.2 

X 

10“ 2 

6.5 

X 

10 _1 

054m 

1.5 

5.550 

X 

10® 

0.75 

1.2 

X 

10“ 2 

7.4 

X 

10“ 2 

4.0 

X 

10 -3 

9.9 

X 

10“ 3 

4.7 

X 

10 _1 

054d 

1.5 

5.550 

X 

10® 

0.75 

1.3 

X 

10“ 2 

7.7 

X 

10“ 2 

2.0 

X 

10“ 3 

1.4 

X 

10“ 2 

8.4 

X 

10- 1 

055d 

1.5 

5.550 

X 

10® 

1.00 

1.3 

X 

10“ 2 

7.7 

X 

10“ 2 

2.1 

X 

10 -3 

1.3 

X 

10“ 2 

7.4 

X 

10 _1 

056d 

1.5 

5.550 

X 

10® 

2.00 

1.2 

X 

10“ 2 

7.4 

X 

10“ 2 




1.5 

X 

1(T 2 

6.3 

X 

10- 1 

057m 

1.5 

6.500 

X 

10® 

0.50 

1.6 

X 

10“ 2 

8.8 

X 

10“ 2 

5.5 

X 

10“ 3 

1.1 

X 

10“ 2 

4.1 

X 

io - 1 

058m 

1.5 

6.500 

X 

10® 

0.75 

1.5 

X 

10“ 2 

8.8 

X 

10“ 2 

4.3 

X 

10 -3 

1.2 

X 

10“ 2 

4.8 

X 

10 _1 

058d 

1.5 

6.500 

X 

10® 

0.75 

1.5 

X 

10“ 2 

8.9 

X 

10“ 2 

2.4 

X 

10“ 3 

1.6 

X 

1(T 2 

8.3 

X 

io - 1 

059d 

1.5 

6.500 

X 

10® 

1.00 

1.5 

X 

10“ 2 

8.8 

X 

10“ 2 

2.5 

X 

10“ 3 

1.6 

X 

10“ 2 

7.7 

X 

IO -1 

060m 

1.5 

8.000 

X 

10® 

0.75 

1.9 

X 

10“ 2 

1.1 

X 

IO -1 

5.4 

X 

10 -3 

1.5 

X 

10“ 2 

3.8 

X 

IO -1 

060d 

1.5 

8.000 

X 

10® 

0.75 

1.9 

X 

10“ 2 

1.1 

X 

nr 1 

3.3 

X 

10“ 3 

2.0 

X 

1(T 2 

8.2 

X 

IO- 1 

061d 

1.5 

8.000 

X 

10® 

1.00 

1.9 

X 

10“ 2 

1.1 

X 

nr 1 

3.0 

X 

10“ 3 

2.0 

X 

10“ 2 

7.7 

X 

IO- 1 

062m 

1.5 

9.000 

X 

10® 

0.50 

2.2 

X 

10“ 2 

1.2 

X 

IO -1 

6.4 

X 

10 -3 

1.6 

X 

10“ 2 

4.5 

X 

IO -1 

063m 

1.5 

9.000 

X 

10® 

1.00 

2.2 

X 

10“ 2 

1.2 

X 

nr 1 

5.1 

X 

10“ 3 

1.8 

X 

1(T 2 

3.1 

X 

IO- 1 

064m 

1.5 

1.000 

X 

10 7 

0.50 

2.5 

X 

10“ 2 

1.3 

X 

nr 1 

8.0 

X 

10“ 3 

1.8 

X 

10“ 2 

3.0 

X 

IO- 1 

065m 

2.0 

3.000 

X 

10® 

1.00 

4.0 

X 

10 -3 

2.9 

X 

10“ 2 

1.8 

X 

10 -3 

3.8 

X 

10“ 3 

2.2 

X 

IO -1 

066d 

2.0 

3.000 

X 

10® 

2.00 

4.0 

X 

10“ 3 

3.1 

X 

10“ 2 

5.5 

X 

10“ 4 

4.5 

X 

1(T 3 

7.7 

X 

IO- 1 

067m 

2.0 

4.000 

X 

10® 

1.00 

6.8 

X 

10“ 3 

4.6 

X 

10“ 2 

1.7 

X 

10“ 3 

6.0 

X 

10“ 3 

2.9 

X 

IO- 1 

068d 

2.0 

4.000 

X 

10® 

2.00 

6.5 

X 

10 -3 

4.6 

X 

10“ 2 

1.1 

X 

10 -3 

7.2 

X 

10“ 3 

7.2 

X 

IO -1 

069d 

2.0 

4.000 

X 

10® 

3.00 

6.6 

X 

10“ 3 

4.6 

X 

10“ 2 

1.1 

X 

10“ 3 

7.3 

X 

1(T 3 

6.4 

X 

IO- 1 

070m 

2.0 

5.000 

X 

10® 

0.50 

8.3 

X 

10“ 3 

5.5 

X 

10“ 2 

5.4 

X 

10“ 3 

6.4 

X 

10“ 3 

2.3 

X 

IO- 1 

071m 

2.0 

5.000 

X 

10® 

1.00 

9.2 

X 

10 -3 

6.0 

X 

10“ 2 

2.1 

X 

10 -3 

7.5 

X 

10“ 3 

2.9 

X 

IO -1 

072d 

2.0 

5.000 

X 

10® 

1.50 

9.3 

X 

10“ 3 

6.4 

X 

10“ 2 

1.6 

X 

10“ 3 

9.9 

X 

1(T 3 

6.8 

X 

IO- 1 

073d 

2.0 

5.000 

X 

10® 

2.00 

9.1 

X 

10“ 3 

6.4 

X 

10“ 2 

1.7 

X 

10“ 3 

1.0 

X 

10“ 2 

6.3 

X 

IO- 1 

074d 

2.0 

5.000 

X 

10® 

3.00 

9.0 

X 

10 -3 

6.1 

X 

10“ 2 

1.5 

X 

10 -3 

1.0 

X 

10“ 2 

6.1 

X 

IO -1 

075m 

2.0 

6.000 

X 

10® 

0.50 

1.1 

X 

10“ 2 

7.3 

X 

10“ 2 

4.1 

X 

10“ 3 

8.2 

X 

1(T 3 

2.5 

X 

IO- 1 

076d 

2.0 

6.000 

X 

10® 

2.00 

1.2 

X 

10“ 2 

7.9 

X 

10“ 2 

2.0 

X 

10“ 3 

1.3 

X 

10“ 2 

5.6 

X 

IO- 1 

077m 

2.0 

7.000 

X 

10® 

0.70 

1.5 

X 

10“ 2 

8.7 

X 

10“ 2 

4.0 

X 

10 -3 

1.2 

X 

10“ 2 

1.7 

X 

10 _1 

078m 

2.0 

7.000 

X 

10® 

1.00 

1.4 

X 

10“ 2 

9.0 

X 

10“ 2 

3.0 

X 

10“ 3 

1.2 

X 

1(T 2 

4.4 

X 

IO- 1 

079d 

2.0 

7.000 

X 

10® 

1.50 

1.4 

X 

10“ 2 

9.0 

X 

10“ 2 

1.9 

X 

10“ 3 

1.6 

X 

10“ 2 

5.9 

X 

IO- 1 

080d 

2.0 

7.000 

X 

10® 

2.00 

1.4 

X 

10“ 2 

8.8 

X 

10“ 2 

1.7 

X 

10 -3 

1.6 

X 

10“ 2 

6.7 

X 

10 _1 

081d 

2.0 

7.000 

X 

10® 

3.00 

1.4 

X 

10“ 2 

9.1 

X 

10“ 2 

1.9 

X 

10“ 3 

1.7 

X 

1(T 2 

5.7 

X 

IO- 1 

082m 

2.0 

8.500 

X 

10® 

0.50 

1.8 

X 

10“ 2 

1.1 

X 

nr 1 

4.6 

X 

10“ 3 

1.3 

X 

10“ 2 

3.9 

X 

IO- 1 

083d 

2.0 

8.500 

X 

10® 

2.00 

1.8 

X 

10“ 2 

1.1 

X 

IO -1 

2.9 

X 

10 -3 

1.9 

X 

10“ 2 

7.5 

X 

10 _1 

084m 

2.0 

1.000 

X 

10 7 

0.50 

2.2 

X 

10“ 2 

1.2 

X 

nr 1 

7.3 

X 

10“ 3 

1.6 

X 

10“ 2 

2.7 

X 

IO- 1 

085m 

2.0 

1.000 

X 

10 7 

3.00 

2.0 

X 

10“ 2 

1.2 

X 

nr 1 

3.9 

X 

10“ 3 

2.1 

X 

10“ 2 

2.1 

X 

IO- 1 

086m 

2.0 

1.200 

X 

10 7 

0.50 

2.6 

X 

10“ 2 

1.4 

X 

IO -1 

8.0 

X 

10 -3 

1.8 

X 

10“ 2 

3.9 

X 

10 _1 

087m 

2.0 

1.400 

X 

10 7 

0.50 

3.0 

X 

10“ 2 

1.6 

X 

nr 1 

1.0 

X 

10“ 2 

2.1 

X 

10“ 2 

3.9 

X 

IO- 1 

088d 

2.5 

3.200 

X 

10® 

4.00 

2.6 

X 

10“ 3 

2.0 

X 

10“ 2 

3.2 

X 

10“ 4 

3.6 

X 

10“ 3 

5.2 

X 

IO- 1 

089d 

2.5 

3.400 

X 

10® 

4.00 

3.1 

X 

10 -3 

2.4 

X 

10“ 2 

3.7 

X 

10 -4 

3.8 

X 

10“ 3 

6.0 

X 

IO -1 

090d 

2.5 

4.400 

X 

10® 

3.00 

5.7 

X 

10“ 3 

4.5 

X 

10“ 2 

8.1 

X 

10“ 4 

7.0 

X 

10“ 3 

5.5 

X 

IO- 1 

091d 

2.5 

4.400 

X 

10® 

4.00 

5.4 

X 

10“ 3 

4.1 

X 

10“ 2 

7.5 

X 

10“ 4 

6.8 

X 

10“ 3 

5.2 

X 

IO- 1 

092d 

2.5 

5.400 

X 

10® 

2.00 

8.0 

X 

10 -3 

5.9 

X 

10“ 2 

1.1 

X 

10 -3 

9.3 

X 

10“ 3 

4.0 

X 

10 _1 

093d 

2.5 

5.400 

X 

10® 

3.00 

7.9 

X 

10 -3 

5.9 

X 

10“ 2 

1.1 

X 

10“ 3 

1.0 

X 

10“ 2 

5.6 

X 

IO- 1 

094d 

2.5 

5.400 

X 

10® 

4.00 

7.6 

X 

10 -3 

5.6 

X 

10“ 2 

1.1 

X 

10 -3 

1.1 

X 

10“ 2 

5.0 

X 

IO -1 

095m 

2.5 

6.400 

X 

10® 

1.00 

9.8 

X 

10 -3 

7.2 

X 

10“ 2 

2.0 

X 

10 -3 

9.0 

X 

10“ 3 

1.3 

X 

10 _1 

096d 

2.5 

6.400 

X 

10® 

2.00 

1.0 

X 

10“ 2 

7.1 

X 

10“ 2 

1.2 

X 

10“ 3 

1.2 

X 

10“ 2 

4.3 

X 

IO- 1 

097d 

2.5 

6.400 

X 

10® 

3.00 

9.7 

X 

10 -3 

6.9 

X 

10“ 2 

1.3 

X 

10 -3 

1.4 

X 

10“ 2 

5.5 

X 

10 _1 

098m 

2.5 

7.400 

X 

10® 

1.00 

1.2 

X 

10“ 2 

8.3 

X 

10“ 2 

2.4 

X 

10 -3 

1.1 

X 

10“ 2 

1.5 

X 

10 _1 

099d 

2.5 

7.400 

X 

10® 

2.00 

1.3 

X 

10“ 2 

8.4 

X 

10“ 2 

1.8 

X 

10“ 3 

1.5 

X 

10“ 2 

6.3 

X 

IO- 1 

lOOd 

2.5 

7.400 

X 

10® 

3.00 

1.2 

X 

10 -2 

8.4 

X 

10“ 2 

1.9 

X 

10 -3 

1.4 

X 

10“ 2 

4.8 

X 

10 _1 

101m 

2.5 

9.000 

X 

10® 

2.00 

1.5 

X 

10“ 2 

9.9 

X 

10“ 2 

3.6 

X 

10 -3 

1.5 

X 

10“ 2 

1.5 

X 

10 _1 

102m 

2.5 

1.000 

X 

10 7 

1.00 

1.8 

X 

10“ 2 

1.3 

X 

nr 1 

4.8 

X 

10“ 3 

1.4 

X 

10“ 2 

1.2 

X 

IO- 1 
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Table Al: (continued) 


model 

N e 

Ra 

Pm 

Ro 

Rot 

Ro z 

Lo 

/diPax 

103m 

2.5 

1.100 x 10 Y 

1.00 

2.4 x 10“* 

1.5 x 10 _1 

5.8 x 10 -3 

1.6 x 10 -2 

1.0 x 10 _1 

104d 

3.0 

8.000 x 10® 

4.00 

9.2 x 10 -3 

6.9 x 10 -2 

1.4 x 10 -3 

1.4 x 10 -2 

6.0 x 10 _1 

105m 

3.0 

9.000 x 10 6 

4.00 

1.2 x 10“ 2 

8.6 x 10“ 2 

3.0 x 10“ 3 

1.4 x 10“ 2 

3.2 x 10 1 
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